Tight lower bound to the geometric measure of quantum discord 
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Dakic, Vedral and Brukner [Physical Review Letters 105,190502 (2010)] gave a geometric measure 
of quantum discord in a bipartite quantum state as the distance of the state from the closest classical 
quantum (or zero discord) state and derived an explicit formula for a two qubit state. Further, 
S.Luo and S.Fu [Physical Review A 82, 034302 (2010)] obtained a generic form of this geometric 
measure for a general bipartite state and established a lower bound. In this brief report we obtain 
' a rigorous lower bound to the geometric measure of quantum discord in a general bipartite state 

which dominates that obtained by S.Luo and S.Fu. 
PACS numbers: 03.67.-a, 03.65.Ta 

Understanding quantum correlations in a multipartite quantum state is a fundamental open prob- 
»—5 1 lem. Over the last two decades quantum correlations are studied in the entanglement-separability 
scenario [l[ leading to important insights in quantum computing quantum communication pro- 
tocols like teleportation [1,01, superdense coding [5], cryptography p etc. However, recently it was 
i — i shown that even some separable states contain nonclassical correlation and can be used to accom- 
plish information processing tasks which cannot be achieved classically_[7-13]. These nonclassical 
correlations of bipartite states are measured by quantum discord 14j, [l5|-the discrepancy between 



quantum versions of two classically equivalent expressions for mutual information. Quantum dis- 
cord became a subject of intense research in different contexts [16-44]. As the evaluation of quantum 
q-i discord involves optimization procedure, analytical results are known only in a few families of two- 



qubit states [18|, [29|, |44|. Recently, a necessary and sufficient condition for the existence of non-zero 



CO 1 quantum discord was obtained and a geometrical way of quantifying quantum discord was proposed 



[451] . Geometric measure of quantum discord, introduced by Dakic et al 45], is given by 



OV D(p) = mm \\p- x\\ 2 , (1) 

where f2o denotes the set of zero-discord states and \\p — x\\ 2 := tr(p — x) 2 is the square norm in the 
Hilbert-Schmidt space. A state x £ H a ® H b is of zero discord if and only if it is a classical-quantum 



state [14j, LlSj] , which can be represented as 



><: x = Y,Pk\k){k\® Pk , (2) 



where {pk} is a probability distribution, {\k)} is an arbitrary orthonormal basis in H a and p k is a 
set of arbitrary states (density operators) acting on H b . 

The quantum discord of a bipartite state p on a system H a <g> H b with marginals p a and p b can be 
expressed as 

Q(p)=min{/(p)-/(n a (p))}. (3) 
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Here the minimum is over von Neumann measurements (one dimensional orthogonal projectors 
summing up to the identity) IP = {11^} on subsystem a, and 

n a ( P ) = ^(n^/ fe )p(n^/ b ) 

k 

is the resulting state after the measurement. I(p) = S(p a ) + S(p b ) — S(p) is the quantum mutual 
information, S(p) = —trplnp is the von Neumann entropy, and I b is the identity operator on H b . 
Intuitively, quantum discord may thus be interpreted as the minimal loss of correlations (as measured 
by the quantum mutual information) due to measurement. This formulation of quantum discord is 
equivalent to the original definition of quantum discord by Ollivier and Zurek [141 ] . 

In order to obtain the desired lower bound for the quantum discord in an arbitrary bi-partite state, 
we set up the following scenario. Consider a bipartite system H a ® H b with dim H a = m and dim 
H b = n. Let L(H a ) be the space consisting of all linear operators on H a . This is a Hilbert space 
with the Hilbert-Schmidt inner product 

(X\Y) = trX ] Y. 

The Hilbert spaces L(H b ) and L(H a £g> H b ) are defined similarly. Let {Xi : % = 1,2, •• • ,m 2 } and 
{Yj : j = 1, 2, • • • , n 2 } be sets of Hermitian operators which constitute orthonormal bases for L(H a ) 
and L(H b ), respectively. Then 

trXiXi' = 5u', trYjYji = Sjjr. 

{Xi<g)Yj} constitutes an orthonormal (product) basis for L(H a ®H b ) (linear operators on H a ®H b ). 
In particular, any bipartite state p on H a (g) H b can be expanded as 



P 



»v-v, y r (4) 



with Cij = tr(pXi <S>Yj). 

S. Luo and S. Fu introduced the following form of geometric measure of quantum discord 1461 
(Here and throughout this article, the superscript t denotes transpose of vectors or matrices). 

D(p) = triCC 1 ) - max tr(iCC*i'), (5) 

A 

where C = [cij] (Eq.(j3])) is an m 2 x n 2 matrix and the maximum is taken over all m x m 2 -dimensional 
isometric (see below) matrices A = [aki] such that 

a ki = tr(\k)(k\Xi) = (k\Xi\k), fc = l,2,...,m 

2 = 1,2,. ..,m 2 , (6) 

and {\k)} is any orthonormal basis in H a . We can expand the operator \k)(k\ in this basis as 

\k)(k\=^2 aki Xi, k = l,2,...,m. (7) 

i 

A = [aki] is an isometry in the sense that AA f = I a and the row vectors 3k of the matrix A satisfy 

2 

m 

I>« = 1- (8) 



|2 _ 

i=l 
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Further, using their definition, it immediately follows that 

m 

y^qfej = trXj. (9) 



fc=i 



We can represent the density operators acting on a bipartite system H a <g) H b , with dimH a = m 
and dimH = n, as 



— {I m <S> I n + ^ ^Ai 8) 7 n + y~l ® A,- 

+ ^^-A i ®A i ), (10) 



where Aj,z = l,...,m 2 — 1 and Aj,j = l,...,n 2 — 1 are the generators of SU(m) and SU(n) 
respectively, satisfying tr(XiXj) = 2Sij [47|. Notice that x E M m _1 and y e R™ ~ x are the coherence 
vectors of the subsystems A and B, so that they can be determined locally. These are given by 



49J,|5QJ 



m ~ m , ~ 

Xi = —tr(pX, t <g> I n ) = —tr(p A Xi 



Tl ~ Tl ~ 



where = ^ r s(p) and p_e = tr^p) are the reduced density matrices. The correlation matrix 
T = [tij] is given by 

r = N = — [^(pAi® A,)]. 

In this article, we find the lower bound of geometric measure of quantum discord which dominates 
the lower bound in ref. |46| . 



Theorem 1. Let p be a bipartite state defined by Eq. (fTUj) ; then 

2 2 m_1 

(pii 2 +-imi 2 -5>), (ii) 

m z n n 

j'=i 

where 77^ , j = 1,2, •• • ,m 2 — 1 are the eigenvalues of the matrix (xx? + 2231) arranged in non- 
increasing order (counting multiplicity). 
We prove this theorem for arbitrary (finite) m and n. 

In Eq.flSj) giving the quantum discord D(p), the maximum in the second term is taken over the 
m x m 2 isometric matrices A which also satisfy Eq.sflH]) and Eq.flT]). In other words, the row vectors 
of A are required to be the coherent vectors of states forming an orthonormal basis in H a . If we 
ignore this constraint while maximizing tr^ACC 1 A 1 ') and maximize over the isometric matrices A 
defined below via Eq. (I18|I19|) . the resulting maximum will be greater than or equal to the required 
maximum of tr^ACC 1 A 1 } over the matrices A satisfying Eq.s( l6|7|) . Since all the terms are positive, 
this leads to 

D{p) > tr(CC l ) - maxtriACC'A 1 ). (12) 
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We proceed to obtain the maximum in the above inequality (Eq. (ll2p ). We choose the ortho normal 
bases {Xi} and {Yj} in Eq.fll]) as the generators of SU(m) and SU(n) respectively (47}. 



X\ — —t=Imi Yl — —J=In 

>m Jn 



and 



Xi = i = 2,3, . . . ,m 2 

V2 



Yj = -s^i-ij 3 = 2,3, . . . ,n 2 . 



y/2 

Since trA; = 0; % — 1, 2, ■ • • , m 2 — 1, we have, via Eq.flH 



a/ci = trX; = trAj = 0,2 = 2,..., m 2 

fe=i 



Therefore, 



m— 1 



a™ = - ^a fci , z = 2,3,- •• ,m 2 . (13) 
fc=i 

We now proceed to construct the m x m 2 matrix ^4 defined via Eq.flBJEJ). We will use Eq. flTBl . The 
row vectors of A are 

0>k = (Qfelj «fc2i ' ' ' , a km 2 )'i k — 1,2, ... ,m. 

Next we define 



4 = J -{a k 2,a k 3, ■ ■ ■ ,<w0> k = 1, 2, . . . , m - 1 (14) 

771 — 1 



and using Eq. flTB"]) . we get 



m—l 



e m = -^2e k . (15) 

k=l 

We can prove 

||e fc || 2 = l k = l,2,...,m-l (16) 

2 

using the condition \\a k \\ 2 = Y^iLi a ki = 1 (Eq.©) and using Eq.flU]) with i — 1, namely, a^i = 
tr(|/c)(A;|Xi) = Further, isometry of the A matrix (AA l = I) implies 

e i e) = ^—,j^i = l,2,...,m-l. (17) 
J m — l 

We can now construct the row vectors of m x m 2 matrix A, using Eq. (jl4j) and Eq. fll 5ft . 

1 



(1, y/m — lik), k = 1, 2, • • • , m — 1 
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m— 1 



^-V^VeJ (19) 



defining matrix A 

We get the elements of C = [cij] = [tr(pXi <g> Yj)] using the definitions of the bases {Xi} and {Y}} 
given above, in terms of the generators of SU(m) and SU(n). This gives 

C = 

and 

fr.(cc*) = (— + -l-IMi 2 + 

mn n m rrrn 

|T|| 2 ). (20) 




n 2 m? 



where x, y and T are the coherent vectors and the correlation matrix respectively, defined in Eq.flTUJ). 
Having constructed the matrices A and C, we get, for tr(ACC t A t ), 



tr(ACC l A u 



11 2 INI , 2fm- 1^ 
m \ n n z 



m 2 n 



m—l m—l m—l 

E ^ + E E 

.i=i *=i i=i 



where 



G = xx 1 + 



ITT 1 



n 



(21) 



(22) 



is the (m 2 — 1) x (m 2 — 1) real symmetric matrix. The eigenvectors of G span IR m 1 and form a 



orthonormal basis of M m . Let 771, 7/2, ... , f]m 2 -i be the eigenvalues of G arranged in non-increasing 
order (counting multiplicity). Let I/2), • • • , \fm 2 -i)) be the corresponding orthonormal eigen- 

vectors of G. To maximize the right hand side of Eq. (l2"T]) . we choose 



ei = fi 



and expand {e} 1 JL 2 in the eigenbasis of G. Thus, 



m—l 



ej = E £ k ] fk j = 2, 3, . . . , m - 1. 



(23) 



(24) 



fc=i 



Eq.f lTB"]) with $i replaced by t\ gives us 



S3) 



, j = 2,3, . . 1. 



m—l 

Since e x = /1 (Eq.(|25|)) and ||e,-|| 2 = 1 for all j (Eq.fllEJ) we get 

m 2 -l 



(25) 



E e fc 

fc=i 



(i) 2 



1. 



(26) 



We substitute Eq.( jHJ) for e 3 - in the expression for tr(ACC*A t ) (Eq. fl2Tj) ) and use Eq.fl25l) to get 
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tr(ACC l A^ 



11 

m n 



n- 



2{m- 1) 

m 2 n 



m 



m — 1 



m—l m — 1 

2 E E 

j=2 k=2 



MY 



Vk 



m—2 m—l m—l 



i=2 j>i k=2 



(27) 



We have to choose vectors j — 1, . . . ,m — 1, that is, the expansion coefficients {e^}; j 



1, . . . , m — 1 



k = l,...,m 2 — 1 in Eg. (I24j) . consistent with Eg. s ( 12511261) . so as to maximize 



tr(ACCM*). First we note that we must use at least (m — 1) eigenvectors fk in Eq. fl24|) to ex- 
pand all of {ej} j = 1, . . . , m — 1, because otherwise {e,,} becomes a set of m — 1 linearly dependent 
vectors in a subspace of dimension less than m — l, in which case the row vectors of matrix A cease 

(2) * 

to be mutually orthogonal. Thus, for m = 3 a choice like £2 = t\ fi yields d\ ■ 0,2 7^ 0. With eigen- 
values {%} k — 1, . . . , m 2 — 1 arranged in non- increasing order, we now see that, for every choice of 

{e^}; j — 1, . . . , m — 1; = 1, . . . , m 2 — 1 we get the following upper bound on tr(ACC"v4*) 



tr^ACC'A 1 ) < 



11 2 2(m-l) 

-\- + —y + 2 

m \ n n z m z n 



m—l / j— 1 m—l 

^».+ 2 E E^ + %E^ 2 



m 

m 2 -l 



j=2 \fc=2 



m—2 m—l / i—1 

+ 2 EE EW»+*E^ 

i=2 jr>i \fc=2 



k 



k=i 



The maximum value of RHS is then obtained by choosing = for all k > j, that is 



(2f 



J>J°A j = 2,3,...,ro-l, 



fc=i 



which leads, using the fact that el = if k > j, to 



(29) 



11 2 2(m-l) 

^ -{- + — \\y\\ + 2 



max[tr(ACC*A t )] 



m—l 



k=2 



m—l 

E4 ; 

,j=k 



m—2 m—l 

+EE'" UJ) 

i=k j>i 



C k e k 



(30) 



The remaining task is to find the coefficients ejj , , fc = 1,2, ... , j in the expansion of ej in the 



eigenbasis {fk}, k = 1,2, . . . , j which satisfy Eq.s fll5|16|r29|) and the consequential equations (I25II26P 
Using Eq.sf ll6|riT|25|26|r29|) the problem can be reduced to the coupled pair of equations 



2 m(m — 2) 



3-1 



(m - iy 



(31) 



k=2 



S3) 



1 

M 



-m 



[m 



/ J t k t k 

k=2 



;j>i> 1. 



(32) 



which can be solved iteratively, starting from j = 2. The result is 

m — j 



ma = to 

3 (m - r 

and 



m — j + 1 



(33) 



A = , =. d>i (34) 

y (to — 1)(to — i + 1)(to — i) 



Using Eq.s (l33|l34l) the last term in Eq. (j3"0]) can be evaluated. We have, 



m— 1 m— 2 m—1 

2(m- 1 



j=k 

Eq.()35D simplifies Eq.flSD to 



E e * + E E W = ^TTTT' < 35 > 

j'=fc i=fc J>i 



max^CC^')] = — + -|-|M| 2 + E>. (36) 
A mn n z m m z n z — ' 

3=1 

Finally, Eq.p]}, Eq.flM]) and Eq.(H2J together imply 

m—1 



i=i 



which completes the proof of the theorem. 
Now, we prove that 



„ „ m—1 m 

_ ... . . r) 



£>(p) > — ai*ii 2 + -imi 2 - E $ ^ tr (^) - E A * ( 37 ) 

m z n n z — ' J z — ' 

j=l i=l 



where Xj are the eigenvalues of CC l listed in decreasing order (counting multiplicity) and r/j are 
the eigenvalues of G = xx 1 + listed in decreasing order (counting multiplicity). The last lower 
bound in Eq. (137j) is proved in ref . |46| where the inequality max^ triACC'A 1 ) < YZi A l is derived. 
Since we have proved max A tr(ACC*A*) < ^ + J^\\y\\ 2 + ^ Vi = a + J^E™iV> to 

prove Eq. (l3"7|) . it is enough to prove 



2 



E A t>«+^E^- (® 

Let us rewrite CC* as 



1=1 8=1 



a it 



M -4-G 



where a = — + \y\ I 2 , G = xx? + 2 -^- and u = -^h=x H — ^f^=Ty. Let eigenvalues of be 
{r^} = ;^{%}) where {7^} are the eigenvalues of G. Then 



2 '"^ 



rrrn 

j=i i=i 



Let \\ < \l < • • • < A a an d Vi — V2 — ■ ■ ■ — 7 2_i- F rom theorem (4.3.8) in [5 lj . we have 

X\ < f]'i < \l < f]2 < ■ ■ ■ < V^2-l — Kn 2 > 

so that 



E >i s E »?• 



(40) 



i=i 



Now, we use Eq. ([3"9~j) to get 



m — 1 



i=m?— m+1 J— 1 j=m 2 — m+1 



(41) 



8=1 



Replacing the first term on the left hand side of Eq. (l4"Tj) by Y^=\ ™ ^7 we using Eq.(l4"0 



m — 1 



E A ?- E 

i=m?— m+1 j=m 2 — m+1 



(42) 



"t" I 

Finally, we use A] = A^ 2 - +1 and 



t 4 

^ = n 



—j 



52j], to get 

m— 1 



i=i 



Examples 

(1) We consider the two qutrit state 



P = P\ e )( e \ + (1 — P) 



9 



(43) 



where \e) = ^(|2) <g> |2) + |3) <g> |3) + |2) ® |1) + |1) ® |2) + |1) ® |3) + |3) ® |1)), I is the identity 
operator and i — 1,2, 3} is the standard basis in C 3 . Fig. 1 shows the variation of lower bound 
(given in Eq.f TTTj) ) and the lower bound on D(p), as given in j46], namely, tr(CC l ) — Y^iLi A* (where 
{Aj} are the eigenvalues of CC* listed in the decreasing order, counting multiplicity), with p. 




Fig. 1 Lower bound of quantum discord (Eq.JTTJ), solid line) and lower bound (Eq.(6) in [4g|, dashed line) as a function of ; 

(Eq.Q3)). 
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We see that the lower bound in Eq. fflTj) dominates this lower bound for p > 0.2. 
(2) We consider the two qutrit state 



P = p|ei)(ei| + (l-p)|e 2 )(e 2 | (44) 

where \ex) = ±|11) + §|22) + ^=|33) and |e 2 ) = ^-(|2) g> |2) + |3) g> |3) + |2) ® |1) + |1) ® |2) + 
|1) <8> 1 3) + 1 3) <8> 1 1) ) . Fig. 2 shows the variation of D(p) (given in Eq. lflT]) and the lower bound on 
D(p), as given in |46|, namely, tr(CC*) — Y1T=\ ^« (where {A,} are the eigenvalues of CC* listed in 
the decreasing order, counting multiplicity), with p. 




0.2 0.4 0.6 0.8 1 

P 

Fig. 2 Quantum discord (Eq.{TT}, solid line) and lower bound (Eq.(6) in _46}, dashed line) as a function of p (Ea. (|44[) 



We see that the lower bound in Eq. lflT]) dominates this lower bound. 

To the best of our knowledge, only the lower bounds on the quantum discord in an arbitrary 
bipartite state are presently available. This seems quite surprising as detecting quantum discord 
is not a NP-hard problem such as the separability problem and one expects an exact computable 
expression for quantum discord in all dimensions. This expectation is further augmented by the fact 
that the set of zero discord states is of measure zero. In order to get such an exact computable 
expression for quantum discord in all dimensions, we have to maximize the second term in Eq.Q 
over the set of m x m 2 matrices A which satisfy (in addition to other conditions) the requirement 
that each row vector of A must be a coherent vector of a orthonormal basis state in H a . For the 
two qubit case, this requirement becomes redundant as every unit vector in K 3 is a coherent vector 
of some single qubit pure state. If dim(H a ) > 2, not every unit vector in R d2 ~ l is a coherent vector 
of some pure state in H a . Hence whenever dim(H a ) > 2, the above requirement is to be included 
as an independent constraint in the constrained optimization of the second term in Eq.flS]). The 
resulting constrined optimization problem is very difficult because the set of coherent vectors (for 
dim(H a ) > 2) do not have some simple geometric structure like Bloch sphere. This is the reason why 
an exact computable expression for quantum discord for all dimensions still eludes us. We note that 
the lower bound in Eq.t fiT]) becomes exact for a2xd system (with measurement on the qubit) since 
the above constraint is relaxed in this case. Finally, our lower bound on discord can be meaningfully 
used to compare discordant states and check on the possible monogamy property of dicord [531 ] . 
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